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ABSTRACT 

We report on new measurements of breakup temperatures for target spectators 
from ^^^Au + ^^'^Au reactions at 1000 MeV per nucleon. The temperatures rise 
with decreasing impact parameter from 4 MeV for peripheral to about 10 MeV for 
the most central collisions, in good agreement with previous results for projectile 
spectators at 600 MeV per nucleon. 

The measured temperatures agree quantitatively with the breakup temperatures 
predicted by the statistical multifragmentation model. For these calculations a re- 
lation between the initial excitation energy and mass was derived which gives good 
simultaneous agreement for the fragment charge correlations. 
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The energy spectra of light charged particles and the behaviour of the mean kinetic 
energies of neutrons indicate a substantial component of light particle emission prior 
to the final breakup stage. 



1 Introduction 



The correlated measurements of the nuclear temperature and the excitation energy 
for excited projectile spectators in ^^^Au + ^^^Au collisions at 600 MeV per nucleon 
has permitted to extend the caloric curve of nuclei to very high excitation energies, 
far beyond the nuclear binding energy 0]. For the temperature determination 
the method suggested by Albergo et al. has been used which is based on the as- 
sumption of chemical equilibrium and requires the measurement of double ratios of 
isotopic yields 0. The obtained temperatures, plotted against the measured excita- 
tion energy, resulted in a caloric curve with the characteristic behavior reminiscent 
of first-order phase transitions in macroscopic systems. The 'liquid' and the 'vapor' 
regimes where the temperature rises with increasing excitation energy are separated 
by a region of nearly constant temperature T ^ 5 MeV over which the multiplicity 
of the fragmentation products increases. 

The present discussion of the caloric curve mainly centers around three points. 
The first one concerns the methodical question of whether breakup temperatures 
can be reliably deduced from isotopic yield ratios in the presence of sequential-decay 
processes by which these ratios may be modified after the breakup has occurred 
0, H. Model calculations show that the amount of modification depends on the 
breakup density which is not very well known [2,6-8]. The potential directions of 
interpretation, in a wider sense, of these new experimental results constitute the 
second topic of discussion [9-11]. Finally, there is numerous work concerned with the 
consequences for finite nuclei [O, O] of the predicted first-order phase transition in 



nuclear matter [0, ^ (for very recent references see, e.g., refs. fT^, [T3])- Iii addition, 
experimental caloric curves for other systems have been reported recently [jl8|, |19|. 

In this contribution, we present results of new measurements of the breakup tem- 
perature in the ^^^Au + ^^^Au reaction. In particular, we will address two specific 
questions connected with the statistical interpretation of multi-fragment decays of 
excited spectator systems: Do the measured temperatures exhibit the same invari- 
ance with respect to the entrance channel that has been found for the fragmentation 
patterns? This universal feature of the spectator decay has become evident in the 



observed Zhound scaling of the fragment multiplicities and charge correlations [20 



Here Zhound, defined as the sum of the atomic numbers Zi of all spectator fragments 
with > 2, is a quantity closely correlated with the excitation energy transferred 
during the initial stages of the reaction. Secondly, how do the extracted tempera- 



tures compare with the predictions of the statistical multifragmentation model [13] 
which has been so successfuU in describing the multiplicities and charge correlations 
characterizing the universal partition space ||21||? 



2 Experimental method and results 



The experiment was performed at the ALADIN spectrometer |2y] of the GSI facihty 
with beams of ^^^Au of 1000 MeV per nucleon incident energy, provided by the heavy- 
ion synchrotron SIS. The present data were taken as part of a larger experiment which 



incorporated three multi-detector hodoscopes for coincident particle detection p, ^ 

A set of seven telescopes, each consisting of three Si detectors with thickness 50, 
300, and 1000 fim and of a 4-cm long CsI(Tl) scintillator with photodiode readout, 
were used to measure the isotopic yields of light charged particles and fragments. 
Four telescopes were placed in the forward hemisphere while three telescopes were 
placed at 9iab = 110°, 130°, and 150° for detecting the products of the target-spectator 
decay. Isotopes in the range from hydrogen to carbon were satisfactorily resolved. 
Because of the presumably different shapes of ^He and ^He spectra at low energies, a 
correction was required in order to compensate for the effect of the detection threshold 
E ^ 12 MeV for helium ions, resulting from triggering with the 300-/im detector. For 
a global characterization of the reaction and impact parameter selection, the quantity 
Abound of the coincident projectile decay was measured with the time-of-flight (TOF) 
wall of the ALADIN spectrometer. 

Emission temperatures T were deduced from the double ratios R formed by com- 
bining the ratio of ^He/^He yields with either the lithium yield ratio ^Li/^Li or with 
the hydrogen yield ratios p/d or d/t. The set of ^He, *^He, ^Li, and ^Li isotopes is 
the one used previously for the determination of breakup temperatures of projectile 
spectators in ^^^Au + ^^^Au at 600 MeV per nucleon |jl], ^j. Combinations involving 
p, d, or t, together with ^He and "^He, have the advantage of larger production cross 
sections, particularly in the 'vapor' regime where the heavier fragments are becoming 
rare. The three expressions are 

THeLi,o = 13.3/ln(2.2^i4^) 

and 

THepd,o = 18.4/ln(5.6^^^) (2) 

y mc/ 4He 

and 

THedt,o = 14.3/ln(1.6^^^i^) (3) 

where the temperatures are given in units of MeV. 

The subscript is meant to indicate that these apparent temperatures, derived 
from the measured ground-state populations, may be affected by feeding of these 
populations from sequentially decaying excited states. The required corrections were 
calculated with the quantum statistical model which starts from chemical equilib- 
rium at a given temperature, density, and neutron-to-proton (N/Z) ratio and which 
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Figure 1: Apparent temperatures TneLi.O) ^Hepd.Oi and THedt.O) according to the quantum 
statistical model, a,s ci function of the input tGnipcrciturc Tifiput- 

A breakup density p/ pQ = 
0.3 is assumed. The dotted line represents the linear relation Tq = Tj„pu(/1.2. 



includes sequential decay p3|, |2J]. In Fig. 1, the three apparent temperatures defined 
in eqs. (1-3) are shown as a function of the equilibrium temperature Tinput for the 
parameters N/Z = 1.49 (value of ^^^Au) and density p = 0.3 ■ po (where po is the 
saturation density of nuclei). The relations between TneLi.o or Tnedt.o and T^put are 
almost linear and the corrections required in these two cases are practically identi- 
cal, except at the highest temperatures. The linear approximation, indicated by the 
dotted line, corresponds to the constant correction factor T = 1.2 ■ Tq adopted for 
TneLi in Ref. [Q. Apparently, Tnepd.o is more strongly affected by feeding effects at 
the higher temperatures. Within the range of densities 0.1 < p/po < 0.5, the correc- 
tions required according to the quantum statistical model vary within about ±15% 
0, . They virtually do not change with the N/ Z ratio of the primary source. In 
the analysis, the corrections calculated for p/po = 0.3, as derived from the results 
shown in Fig. 1, were applied. 

In Fig. 2 the obtained temperatures TneLi, ^Hepd, and Tnedt are shown as a 
function of Zbound- They are based on the yield ratios measured with the telescope 
at the most backward angle 6iab = 150°. Simulations indicate that, at this angle, 
contributions from the midrapidity source are negligible. The temperatures increase 
continuously with decreasing Zbound from T = 4 MeV for peripheral collisions to 
about 10 MeV for the most central collisions associated with the smallest Zbound 
values. The range Zbound < 20 corresponds to the high excitation energies at which 
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Figure 2: Temperatures TkeLii ^Hepd, and Tnedt as a function of Zhound, averaged over 
bins of 10-units width. Corrections have been apphed as described in the text. The 
error bars represent the statistical uncertainty. The systematic uncertainty, caused by the 
extrapolation of the yields of helium isotopes below the identification threshold, is indicated 
by the brackets. For clarity, the open data symbols are laterally displaced by 2 units of 

Zbound ■ 



the upbend of the temperature appears in the caloric curve |1| . The results obtained 
with the three different double ratios agree rather well which is remarkable in view of 
the different feeding corrections that are required, in particular, in the case of Tnepd- 



3 Universality of spectator temperatures 

The universality of spectator fragmentation has first been recognized in the study 



of the charge partitions [^, |26[. Universality, in the present context, refers to the 
invariance of the fragmentation patterns of excited spectator nuclei with respect to 
the properties of the entrance channel. The decaying system has lost the memory of 
how it was formed. This was shown to be valid for the variation of the bombarding 
energy (beyond about 400 MeV per nucleon) and for the variation of the target 
mass (if we consider the case of projectile fragmentation). If a linear scaling law is 
applied, it is even valid for the variation of the projectile mass itself. The search 
for the deeper symmetries behind these universal properties of the fragment decay is 



presently pursued by several groups p^, EH 



The data show that the invariance with respect to the entrance channel includes 




Figure 3: Temperatures ThgLI of the target spectator from the present experiment at E/A 
= 1000 MeV (dots) and of the projectile spectator (from Ref. |2§]) at E/A = 600 MeV 
(open squares) and 1000 MeV (open circles) as a function of Zf^gy^^^. The data symbols 
represent averages over bins of 10-units width. For clarity, the open data symbols are 
laterally displaced by 2 units of Z hound- The full and dashed lines represent the breakup 
temperature Tbreakup and the isotopic temperature TneLi calculated with the statistical 
multifragmentation model. 



the deduced isotope temperatures. In Fig. 3 the TneLi temperatures from this work 
and those derived previously [|l], |29| for projectile spectators in ^^''Au + ^^^Au col- 
lisions at 600 and 1000 MeV per nucleon are shown in comparison. In the case of the 
projectile decays, the isotopes were identified by tracking of their trajectories with the 
upgraded TP-MUSIC detector and subsequent momentum and time-of-flight analy- 
sis. The agreement obtained in the measurements at 1000 MeV per nucleon for the 
projectile and the target decays, performed with different methods of isotope identi- 
fication and associated with different detection thresholds, reflects the experimental 
accuracy. The agreement of the data measured at 600 and 1000 MeV per nucleon 
confirms the expectation that the breakup temperature, as a function of Zhound, does 
not change with the bombarding energy. 

Neither does the breakup temperature depend on the mass of the collision partner. 
In Fig. 4 this is demonstrated for the case of ^^^Au projectiles of 600 MeV per nucleon 
impinging on C, Al, Cu, Au targets. The breakup temperatures deduced by the EOS 



collaboration for Au + C at 1 GeV per nucleon |18] are also consistent with this 



conclusion. Within the range Zhound > 40 which is mainly populated in this reaction 
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Figure 4: Top: Measured cross sections da / dZ hound for the reactions of ^^'''Au projectiles 
at E/A = 600 MeV with the four targets C, Al, Cu, and Au. Note that the experimental 
trigger affected the cross sections for ^feound 

> 65. 

Bottom: Breakup temperature TueU for the same four reactions (from Ref. ||25[| ). 



(Fig. 4, top, and Ref. |]20|), they are in good agreement with the present results 
for ^^^Au + ^^^Au, both in absolute magnitude and in their dependence on the 
impact parameter. The observed universality is a strong indication for equilibrium 
being reached at the breakup stage which is a necessary condition for a statistical 
description of the fragmentation process. 

4 Statistical model description 



Calculations within the statistical multifragmentation model [|I3| were performed in 
order to test its consistency with respect to the statistical parameters and predicted 
charge partitions. In the model one assumes that all observed particles come from 
the decay of one equilibrated source. 

As a criterion for selecting the distribution in excitation energy and mass of the 
ensemble of excited spectator nuclei, required as input for the calculations, we chose 
the capability of the model to simultaneously describe the correlations of the mean 
multiplicity (Mjmf) of intermediate-mass fragments (IMF's) and of the mean charge 
asymmetry (012) with Z^ound- The asymmetry an of the two largest fragments is 
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Figure 5: Mean multiplicity of intermediate-mass fragments (Mjmf) (top) and mean 
charge asymmetry (012) (bottom) as a function of Z\,oundi as obtained from the calcula- 
tions with the statistical multifragmentation model (open circles) in comparison to the 
experiment (dots, from Ref. |^]). The dashed and dotted lines show the results of the 
calculations with excitation energies E^/A 15% above and 15% below the adopted values, 
respectively. Note that the trigger threshold affected the data of Ref. ||2^ at Zhound > 65. 



defined as = {Z^ax — Z2)/{Zmax + Z2), with the mean value to be calculated 
from all events with Z^ax ^ Z2 > 2. The comparison, shown in Fig. 5, was based 
on the data reported in Ref. [^]. In the region Zbound > 30 the mean excitation 
energy of the ensemble of spectator nuclei was found to be well constrained by the 
mean fragment multiplicity alone. At Zbound ~ 30 and below, the charge asymmetry 



was a necessary second constraint (cf. Ref. |]SD[) while, at the lowest values of 
Zbound, neither the multiplicity nor the asymmetry provided rigid constraints on the 
excitation energy. These sensitivities are illustrated in Fig. 5 where the dashed 
and dotted lines show the model results for E^/A chosen 15% above and below the 
adopted values. 

The solid line in Fig. 3 represents the thermodynamical temperature T obtained 
from the calculations. Within the given experimental and methodical uncertainties, 
it is in very good agreement with the data. The description of the fragmentation as a 



statistical process is thus internally consistent, the temperatures needed to reproduce 
the observed partition patterns correspond to those measured. The comparison is 
even more consistent if we take directly the temperature TneLi that is obtained from 
the calculated isotope yields and corrected in the same way as the experimental data 
(dashed line). The agreement of the measured and model Thsli is excellent, and 
the tendency of being lower (higher) than the thermodynamical temperature at the 
larger (smaller) values of Zhound is common to both. According to the statistical 
multifragmentation model, therefore, the true source temperature varies only from 
T = 5 MeV for peripheral to about 9 MeV for the most central collisions and stays 
rather constant T ^ 6 MeV over a wide range of Zhound, i-e. over a wide range of 
excitation energies, where the maximum production of intermediate-mass fragments 
is observed. 

The excitation energies that have resulted from this procedure are somewhat 
larger than what was found previously in analyses [21,30-32] of the earlier ^^''Au on 
Cu data at 600 MeV per nucleon ||33|. The difference reflects the sensitivity to the 



fragment multiplicity and is caused by the slightly larger mean multiplicities that were 
obtained from the more recent experiments with improved acceptance pO|. They are 
still considerably lower than the energies obtained from the calorimetric measurement 
of the total energy deposit (see section 6). A (partial) resolution of this discrepancy 
may come from pre-breakup emission of light particles carrying away energy as the 
spectator system approaches the final breakup stage. 



5 Indications of pre-breakup emission 

Kinetic energy spectra were studied for light charged particles up to ^He. For the five 
species proton, deuteron, triton, ^He, and ^He, the spectra measured at 9iab = 150° 
and integrated over finite regions of Zhound are shown in Fig. 6. For comparison, 
the predictions calculated with the statistical fragmentation model are also given. 
The two sets of experimental and model spectra are each normalized, and one overall 
normalization factor is used to relate the two sets. It was adapted to the yields of Z = 
2 fragments because their calculated multiplicities, as a function of Zhound, are found 
to satisfactorily reproduce the experimental multiplicities reported in Ref. . 



The main trend apparent from the comparison is a systematically increasing de- 
viation of the experimental from the model spectra with decreasing Zhound, i-e- in- 
creasing centrality, and with decreasing particle mass. It not only affects the slope 
parameters describing the shape of the spectra but also the integrated intensities. 
The yields of hydrogen isotopes, and in particular of the protons, are grossly under- 
estimated by the statistical multifragmentation model. In the case of ^He, on the 
other hand, the equilibrium description accounts rather well for the multiplicities 
and kinetic energies. A major contribution to the observed ^He yields is expected 
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Figure 6: Energy spectra, measured at 6iab = 150°, of light charged particles p, d, t, ^He, 
and ^He for three intervals of Zbound as indicated. The dots represent the measured spectra, 
the open circles are the results of the calculations with the statistical multifragmentation 
model. The spectra are normalized as stated in the text. 



to come from evaporation by large fragments and excited residue-like nuclei which, 
apparently, is modelled well. 

Conceivable mechanisms that cannot explain the observed deviations include col- 
lective flow and Coulomb effects which both should act in proportion to the mass 
or charge of the emitted particle, contrary to what is observed. On the other hand, 
the commonly adopted scenario of freeze-out after expansion involves a pre-breakup 
phase during which the system cools not only by adiabatic expansion but also by the 
emission of light particles, predominantly nucleons but also light complex particles 



T0| , |31| , p4| . The spectra should reflect the higher temperatures at the earlier stages 
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Figure 7: Mean kinetic energy of neutrons (full symbols) in the rest frame of the projectile 
spectator as obtained in measurements with LAND for the three bombarding energies 
600 (triangles), 800 (squares), and 1000 MeV per nucleon (circles). For the case of 1000 
MeV per nucleon, a comparison is made between the measured mean kinetic energies of 
protons (open circles) and those obtained by adding estimated Coulomb contributions to 
the neutron kinetic energies (full line). 



of the reaction, prior to the final breakup into fragments, but not necessarily exhibit 



a clear preequilibrium character (cf. Ref. |35 



A significant component of pre-breakup emission in the light particle yields has 
two consequences that deserve particular attention. The pre-breakup yields of pro- 
tons, deuterons, and tritons are included in the double ratios used to determine the 
temperatures Tnepd and Tnedt- This violates the requirement of thermal and chemical 
equilibrium, which is the basic assumption of the method, and thus may shed doubt 
on the meaning of the consistency exhibited by Fig. 2. On the other hand, the quan- 
tum statistical model predicts that, in particular, the p/d ratio varies sufficiently 
slowly with temperature, such that the overall ratio will not be significantly affected 
by contributions from higher temperatures. The second point concerns the interpre- 
tation of the excitation energy that is obtained in a calorimetric measurement. 



6 Excitation energy of primary spectators 

Rather small fractions of the initial bombarding energy are imparted to the spec- 
tator nuclei in relativistic collisions. The actual amount of energy deposition can 
only be reconstructed from the exit-channel configuration which requires a complete 
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Figure 8: Reconstructed average excitation energy (Ex) /{A) of the decaying spectator 
system (circles) and mean excitation energy used in the calculations with the statisistical 
multifragmentation model as a function of ^boMnd- 



knowledge of all decay products, including their atomic numbers, masses, and kinetic 
energies. A rather large fraction of it resides in the kinetic energies of the produced 
nucleons and in their separation energies. 

For the analysis for the ^^^Au + ^^^Au reaction at 600 MeV per nucleon, the data 
on neutron production measured with LAND were available On the other 



hand, the hydrogen isotopes were not detected in this experiment and assumptions 
concerning their kinetic energies as well as the p/d/t intensity ratios and the overall 
N/Z ratio of the spectator had to be made. From the present experiment at 1000 
MeV per nucleon, data on hydrogen emission from the target spectator source have 
been obtained. Together with the neutron data measured with LAND at the same 
bombarding energy in the previous experiment, a complete picture on light particle 
emission has become available for this reaction. It was found, in particular, that by 
adding a Coulomb component to the neutron kinetic energies, an assumption that had 
to be made previously, the proton kinetic energies are considerably underestimated 
(Fig. 7). Furthermore, the full analysis of neutron emission at three bombarding 
energies has revealed a significant dependence of the neutron kinetic energies on the 
bombarding energy. This phenomenon contrasts the universal properties of fragment 
production but is clearly consistent with a pre-breakup component of nucleon 
emission exhibiting some residual memory of the initial stages of the reaction. 

In Fig. 8 the result of the present analysis of the spectator energy and the 
excitation energies used for the calculations with the statistical multifragmentation 



model are shown. In both cases, the values have increased considerably, compared to 
those reported previously ||20|. The reasons have been given already. The masses Aq 
of the primary spectator systems, on the other hand, have remained the same as they 
consist mainly of the sum of the masses of the produced fragments of intermediate 
mass. They are in rather good agreement with the predictions of the geometrical 
participant-spectator model p^. The maximum number of fragments, observed at 



Zbound ~ 40, is associated with initial excitation energies of {E^)/ {A) ^ 12 MeV and, 
according to the model calculations, with excitation energies of (Ex) / (A) ^ 6 MeV 
in the equilibrium stage at breakup. 



7 Summary 

Breakup temperatures TneLi, ^Hepd, and Tnedt were measured for target spectators 
in ^^^Au + ^^^Au collisions at 1000 MeV per nucleon. In these reactions multifrag- 
mentation is the dominant decay channel of the produced spectator systems over a 
wide range of excitation energy and mass. The corrections for sequential feeding of 
the ground-state yields, based on calculations with the quantum statistical model, 
resulted in mutually consistent values for the three temperature observables. 

With decreasing Zhound, the obtained temperatures increase from T = 4 MeV 
for peripheral collisions to about 10 MeV for the most central collisions. Within 
the errors, these values are in good agreement with those measured previously with 
the ALADIN spectrometer for projectile spectators in the same reaction at 600 and 
1000 MeV per nucleon. The agreement of the temperatures measured at 600 and 
1000 MeV per nucleon confirms the expected invariance of the breakup temperature 
with the bombarding energy. It is consistent with the observed Zbound scaling of 
the mean fragment multiplicities and charge correlations and supports the statistical 
interpretation of the multi-fragment decay of highly excited spectator nuclei. 

The comparison with the results of calculations within the statistical multifrag- 
mentation model shows that a good simultaneous agreement for the charge partitions 
and for the breakup temperatures can be achieved. A necessary requirement for a 
consistent statistical description of the spectator fragmentation is thus fulfilled. The 
model results also suggest that the true source temperature varies slightly less, from 
T = 5 MeV for peripheral to only about 9 MeV for the most central collisions, and 
stays rather constant T ^ Q MeV over the range of excitation energies where the 
maximum production of intermediate-mass fragments is observed. 

The systematic behavior of the kinetic-energy spectra of light charged particles 
indicates contributions from light-particle emission prior to the final breakup stage. 
This is supported by the variation with bombarding energy of the kinetic-energy spec- 
tra of neutrons. A more quantitative understanding of the role of the pre-breakup 
processes will be essential for the interpretation of temperatures obtained from light- 



particle yields as well as of the excitation energies obtained from calorimetric mea- 
surements of the spectator source. 
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